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EXECUTIVE  SUMMARY 


The  objective  of  this  effort  was  to  verify  a  newly  developed  closed-form  solution  to  the 
transient  response  of  a  submerged  prolate  spheroidal  shell  structure  under  the  loading  of  an  underwater 
shock. 


The  closed-form  solution  for  the  response  of  a  spherical  shell  had  been  {neviously  developed 
(Ref.  1, 2).  Numerical  results  had  been  generated  fw  a  particular  test  case.  These  results  had  been  verified 
by  an  underwater  shock  tesdng  program  (Ref.  3). 

The  current  research  has  attempted  to  extend  these  results  to  the  prolate  spheroidal  sheU  case. 
A  prolate  spheroidal  shell  with  a  majcM’  diameter  approximately  10  times  the  minor  diameter  is  a  reasonable 
approximation  to  the  slu^  of  a  modem  submarine.  A  closed-form  solution  for  this  case  would  be  useful 
to  cross  check  numerical  or  finite  element  solutions  for  the  underwater  shock  response  of  submarine  hulls. 

Graduate  student  research  sponsored  by  die  The  Charles  Stark  Draper  Laboratory,  Inc.,  herein 
after  known  as  Draper,  under  Independent  Research  and  Development  (IR&D)  funding  had  developed  the 
theoretical  basis  for  the  solution  to  this  problem  (Ref.  4).  The  purpose  ot  this  effort  was  to  obtain  specific 
numerical  solutions  for  comparison  with  other  numerical  methods  or  test  data. 

Numerical  solutions  were  attempted  fw  three  test  cases.  The  fust  case  was  an  almost  perfect 
sphere.  This  case  produced  exactly  the  same  results  as  previously  reported  for  the  sphere,  thus  providing  a 
first  cross  check  of  the  theory. 

The  second  test  case  was  an  egg  shaped  prolate  spheroidal  shell.  This  case  gave  good  results 
for  the  natural  frequencies  of  the  shell  in  vacuum,  but  did  not  give  reasonable  results  of  the  shock  loading. 
The  third  case  was  a  submarine  shaped  prolate  spheroidal  shell.  For  this  case,  all  the  numerical  results  were 
obviously  in  error. 

The  theory  appears  to  be  complete  and  correct.  However,  the  answo-s  are  in  the  form  of 
infinite  series  solutions.  In  some  cases,  the  series  are  cross  coupled,  making  for  an  infmite  square  matrix  of 
terms.  For  practical  solution  with  currently  reasonable  computer  resources,  the  cross-coupling  terms  must 
be  dropped  and  the  series  truncated  at  a  reasonable  number  of  terms.  Unfortunately,  this  results  in  incorrect 
answers.  In  the  future,  increased  computer  power  may  make  this  solution  useful  in  an  engineering  sense. 
However,  at  this  time  the  method  can  not  be  used  to  solve  practical  problems. 
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Since  reasonable  numerical  solutions  could  not  be  obtained  for  high  eccentricity  geometries 
ndng  the  new  closed-form  theory,  we  were  unable  to  compare  the  new  thetny  with  experimental  ot  flnite 
element  results. 


Editorial  Note:  The  above  view  does  not  necessarily  agree 
with  that  of  Dr.  Janet  Jones-Oliveira  who  developed  the  theory 
as  presented  in  Reference  4. 
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CONVERSION  TABLE 


angstrom 

1.000  000  X  E  -10 

meters(m) 

atmosphere 

1.013  25  X  E  +2 

kilo  pascal  (kPa) 

bar 

1.000  000  X  E  -*-2 

kilo  pascal  (kPa) 

barn 

1.000  000  X  E  -28 

meter^  (m^) 

British  thermal  unit 
(thermochemical) 

1.054  350  X  F  +3 

Joule  (J) 

calorie  (thermochemical) 

4.184  000 

Joule  (J) 

cal(thermochemical)/cm^ 

4.184  000  X  E  -2 

mega  joule/m^  (MJ/m^) 

curie 

3.700  000  X  E  ^1 

•giga  becquerel  (GBq) 

degree  (angle) 

1.745  329  X  E  -2 

radian  (rad) 

degree  Fahrenheit 

tg  =  (t*  f  +  459.67)/1.8 

degree  kelvin  (K) 

electron  volt 

1.602  19  X  E  -19 

joule  (J) 

•rg 

1.000  000  X  E  -7 

Joule  (J) 

erg/second 

1.000  000  X  E  -7 

watt  (W) 

foot 

3.048  000  X  E  -1 

meter  (m) 

foot-pound'force 

1.355  818 

joule  (J) 

gallon  (U.S.  liquid) 

3.785  412  X  E  -3 

meter^  (m^) 

inch 

2.540  000  X  E  -2 

meter  (m) 

jerk 

1.000  000  X  E  -^9 

Joule  (J) 

joule/kilogram  (J/kg) 
(radiation  dose  absorbed)- 

1.000  000 

Gray  (Gy) 

kilotons 

4.133 

terajoules 

kip  (1000  Ibf) 

4.448  222  X  E  -•-3 

newton  (N) 

kip/inch^  (ksi) 

6.894  757  X  E  +3 

kilo  pascal  (kPa) 
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CONVERSION  TABLE  (Continued) 


ktap 

newton-second/m^ 

1.000  000  X  E  -1^2 

(N-s/m^) 

micron 

1.000  000  X  E  -6 

meter  (m) 

mil 

2.540  000  X  E  -5 

meter  (m) 

mile  (international) 

1.609  344  X  E  -f3 

meter  (m) 

ounce 

2.834  952  X  E  -2 

kilogram  (kg) 

pound-force 

4.448  222 

newton  (N) 

(lbs  aviordupois) 

pound-force-inch 

1.129  848  X  E  -1 

newton-meter  (N-m) 

pound-force/inch 

1.751  268  X  E  +2 

Newton/meter  (N/m) 

pound-force/foot^ 

4.788  026  X  E  -2 

kilo  pascal  (kPa) 

pound-force/inch^  (psi) 

6.894  757 

kilo  pascal  (KPa) 

pound-mass 

4.535  924  X  E  -1 

kilogram  (kg) 

(Ibm  avoirdupois) 

pound-mass-foot^ 

kilogram-meter^ 

(moment  of  inertia) 

4.214  011  X  E  -2 

(kg  -  m2) 

pound-mass/foot^ 

kilogram/meter^ 

1.601  846  X  E  -i^l 

(kg/m3) 

rad  (radiation  dose 

1.000  000  X  E  -2 

••Gray  (Gy) 

absorbed) 

roentgen 

coutomb/kilogram 

2,579  760  X  E  -4 

(C/kg) 

shake 

1.000  000  X  E  -8 

second  (s) 

slug 

1.459  390  X  E  4l 

kilogram  (kg) 

torr  (mm  Hg,  0*C) 

1J33  22  X  E  -1 

kilo  pascal  (kPa) 

The  becquerel  (Bq)  is  the  SI 

unit  of  radioactivity;  1  Bq 

=  1  event/s. 

**  The  Gray  (Gy)  is  the  SI  unit  of  absorbed  radiation 
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SECTION  1 


INTRODUCTION 


Stnictural  acoustics  is  an  area  of  mathematical  physics  which  addresses  the  covqtled  fluid-solid 
interaction  between  a  structure  and  the  fluid  in  which  it  is  immersed.  In  fluid-solid  interaction  problems, 
the  analysis  of  the  response  of  structures  is  coupled  through  the  boundary  conditions  to  the  propagation  oX 
energy  in  the  fluid  environment.  It  is  incorrect  to  analyze  an  immersed  structure  independently  of  the  fluid 
medium,  to  analyze  the  fluid  independently  of  a  structure  in  its  midst,  and  then  to  superimpose  the 
solutions.  These  problems  require  the  simultaneous  solution  of  the  coupled  field  ^uatitMis  for  both  media 
(Ref.  4,  Chapter  1). 

Because  of  the  coupled  nature  of  the  problem,  very  few  closed  form  solutions  have  been 
developed.  The  problem  of  the  shock  response  for  a  sphere  in  water  was  previously  solved  by  Draper 
(Ref.  1, 2).  Numerical  results  were  develtqted  for  a  spherical  shell  structure  which  DrapCT  was  developing 
ftv  a  different  program.  The  numerical  results  were  then  compared  with  actual  shock  test  results  (Ref.  3). 
The  agreement  was  excellent 

Encouraged  by  these  results,  we  attempted  to  extend  the  analysis  to  a  prolate  q)heroidal  shell 
geometry.  Unfortunately,  this  simple  change  in  the  geometry  caused  a  very  large  increase  in  the  difficulty 
of  the  solution.  Without  the  aid  of  computers,  solution  would  be  practically  impossible. 

However,  computer  programs  are  now  available  to  help  with  the  algebra  involved  with 
manipulating  the  complex  series  expressions  which  result  ffom  attempts  to  solve  this  problem  in  closed 
form.  In  particular,  the  DOE-MACSYMA  program  was  used  extensively  to  help  in  reducing  the  equations 
to  more  reasonable  forms. 
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SECTION  2 


THEORETICAL  DEVELOPMENT 


The  theoretical  devek^ent  of  the  equations  is  given  in  detail  in  Jones-Oliveira  (Ref.  4).  We 
will  present  here  only  a  brief  summary  of  the  steps  involved. 

The  objective  is  to  analyze  the  transient  response  of  a  prolate  spherical  shell  structure  which 
has  been  externally  loaded  by  an  end-on  explosive  watobtxne  shock  wave.  The  intention  is  to  analyze  the 
effects  of  giving  up  one  level  of  geometric  symmetry  (two  degrees  of  symmetry  for  the  prolate  spheroid 
versus  three  degrees  of  symmetry  for  the  previous  ^here  analysis)  (Ref.  4,  Chapter  2). 

Several  simplifying  assumptions  are  necessary  to  formulate  a  solution  to  this  problem.  It  is 
assumed  that  linear  acoustic  wave  theory  is  sufficient  to  describe  the  jvopagation  of  the  pressure  Held  in  the 
water.  This  is  a  reasonable  assumption  except  very  near  the  explosive  charge.  Koiter-Sanders-Budiansky 
theory  for  thin  elastic  shells  is  used  for  the  formul^on  of  the  shell  structure  internal  strain  and  kinetic 
energy  expressions  (Ref.  4,  Chapter  2).  It  is  assumed  that  the  fluid  remains  in  contact  with  the  shell  at  all 
times  (no  cavitation).  This  will  be  true  fw  reasonable  explosive  loads  except  very  near  the  surface. 

With  these  assumptions,  the  problem  is  linear.  One  set  trf  diffoential  equations  describes  the 
motion  of  the  fluid.  A  second  set  of  differential  equadcms  describes  die  motion  of  the  shell  structure.  The 
two  sets  of  equations  are  coupled  by  the  requirement  that  the  normal  velocity  and  the  pressure  load  are 
matched  at  all  times  at  the  fluid-solid  interface.  The  problem  has  thus  been  reduced  to  finding  the  solution 
to  these  equations,  a  far  from  trivial  task. 

The  problem  is  formulated  in  a  prolate  spheroidal  coordinate  system.  The  temporally  and 
^Mitially  dependent  shell  displacement  and  fluid  pressure  fields  are  expressed  modally  in  terms  of  the  prolate 
spheroidal  angular  functions  of  the  first  kind  and  a  new  set  of  predate  spheroidal  radial  functions  of  the  first 
and  third  kinds.  The  prolate  spheroidal  angular  functions  are  shown  to  be  solutions  of  both  the  acoustic 
wave  equation,  which  governs  the  fluid  behavior,  and  the  shell  equilibrium  equations.  They  are  not  the 
eigenfunctions  of  the  coupled  fluid-structure  problem.  However,  the  solutions  are  exact  within  the  limits  of 
using  a  finite  numbers  of  terms  of  a  infinite  series  to  represent  the  solution  (Ref.  4,  Abstract). 

The  prolate  spheroidal  special  functions  are  expressed  in  terms  of  classical  spherical 
harmonics:  Specifically,  Legendre  polynomials  and  modified  spherical  Bessel  functions  of  the  first  and 
third  kinds  (Ref.  4,  Abstract). 


2 


The  mathematical  physics  governing  the  symmetry  order  two  prolate  spheroidal  problem  is 
fundamentally  different  from  that  of  the  symmetry  order  three  spherical  problem  previously  solved.  The 
underlying  geometric  differences  are  at  the  root  of  the  modal  coupling,  which  is  compounded  by  the  fluid- 
solid  interaction  at  the  interface.  The  totally  uncoupled  spherical  response  is  a  limiting  case  of  the  more 
general  and  fully  coupled  prolate  q)heroidal  response  (Ref.  4.  Chapter  2). 

The  modal  coupling,  which  is  not  present  in  the  spherical  case,  is  the  primary  limitation  to 
obtaining  practical  solutions.  If  the  coupling  were  not  present,  it  would  be  reasonable  to  carry  a  large 
number  of  terms  in  the  series  solutions.  If  the  series  converged  very  rapidly,  it  would  not  be  too  difficult  to 
solve  the  fully  coupled  problem  using  only  a  few  terms  in  the  series.  Unfortunately,  it  ^[>pears  that  neither 
is  true. 


DOE-MACSYMA  was  used  extensively  to  determine  the  solution  to  this  problem.  The 
capabilities  of  the  symbolic  computer  code  have  been  demonstrated  to  be  invaluable  for  parametric  analysis 
(Ref.  4,  Abstract).  Unfortunately,  the  complexity  of  the  resulting  coupled  functions  is  overwhelming,  even 
with  DOE-MACSYMA.  It  became  obvious  that  simplifrcations  were  imperative  to  achieve  solutions 
within  the  limitations  of  DOE-MACSYMA  and  reasonable  total  computer  time.  It  was  necessary  to 
estimate  the  relative  magniuides  of  each  of  the  coupling  terms  to  determine  when  terms  could  be  neglected. 
These  simplifications  must  be  considered  judgement  calls. 

The  coupling  terms  are  more  significant  as  the  aspect  ratio  is  increased.  In  {nincipal,  this  is 
not  a  limitation.  As  the  aspect  ration  increased,  one  would  just  carry  more  terms  in  the  solution.  In 
practice,  this  is  subject  to  limitations: 

1.  Numerical  precision  of  the  calculations  must  also  be  increased  to  avoid  round-off 
error. 

2.  Computer  programs  have  size  limits.  In  particular,  DOE-MACSYMA  can  handle 
large,  but  not  unlimited,  problems.  This  analysis  was  stretching  its  limits. 

3.  Computers  have  become  faster  and  less  expensive.  However,  problems  which  are 
attempted  have  also  grown.  This  problem  might  have  been  possible  to  solve  to 
higher  accuracy  by  the  use  of  much  larger  computer  resources,  but  it  would  have 
been  prohibitively  expensive  to  try. 
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SECTION  3 


^  Numerical  calculations  are  ptesenied  for  diree  representative  prolate  ^heroidal  geometries.  The  first 
case  models  a  “nearly  spherical”  geometry  and  is  offered  as  verification  of  the  correctness  of  the  code  by 
comparing  the  results  to  the  sperical  analog.  The  second  case  models  a  “lounded-off  football.”  Preliminary 
results  are  offered  to  a  third  case  which  models  a  “submarine-like”  structure.  The  numerical  complexities 
ate  enormous.  Limitations  and  comments  are  noted. 

^  Formulas  to  the  explosive  characteristics  are  given  in  Jones-Oliveira  (Ref.  4,  Chapter  5.4). 
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Table  1. 


Physical  ffwistanty  This  taUe  provides  the  physical  material  properties  used  to  generate  the 
numerical  examples. 


SYMBOL 

DEFINITION 

VALUE 

exponential  decay  constant 

0.3466  ms 

Cf 

sound  velocity  in  fluid 

1498“- 

s 

Cs 

sound  velocity  in  shell 

5404“- 

s 

C* 

mtdimensianalized  sound  ^>eed  squared 

13.0 

E 

Young's  modulus 

2.11  X  10*®-^ 

V 

Poisson's  ratio 

0.30 

P 

peak  external  pressure  given 

1,883.966 

m* 

60  #  of  HBX-1  at  25  ft  S.OJ>. 

n 

nondimensionalized  peak  pressure 

8.2406x10'^ 

a 

density  of  water 

1.0188  X  10* 

m^ 

A 

density  of  shell 

7.9385  X  10* 

m 

specific  weight  of  water 

999.1-^ 

m^ 

AS 

specific  weight  of  shell 

7785-^ 
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The  geometric  parameters  associated  with  each  of  the  three  cases  are  given  in  Table  2  on  page 
9.  In  each  case,  a  specific  shell  geometry  was  selected  and  sized,  and  the  shell  thickness  was  determined  to 
ensure  a  neutrally  buoyant  body. 

The  neutral  buoyancy  condition  is  defined  such  that  the  mass  of  the  fluid  displaced  by  the 

prolate  ^betoid  must  equal  liie  mass  of  the  prolate  spheroidal  shell.  The  volume  of  a  isolate  spheroid  is 

givoi  by  V  ■  ^aab^  =  ^  rr  |  -l)  (Ref.  4,  Chapter  2.4.1).  Clearly,  the  neutral  buoyancy  condition  is 

3  3 

a  function  of  the  nondimensionalized  mass  ratio,  but  it  is  also  a  function  of  the  eccentricity.  The 
nondimensionalized  shell  thickness  may  be  determined  by  solving  the  following  cubic  equation. 


999.l[(fc*i)’.(s, 


The  remaining  geometric  parameters  are  then  nondimensionalized  as  per  (Ref.  4,  Chapter  2.7). 
Note  that  the  nondimensionalizations  of  time  differ  for  each  of  the  three  cases;  therefore,  so  must  the  scaled 
values  of  the  incident  pressure  exponential  decay  constant 

Case  I:  Model  of  a  nearly  spherical  shell.  The  prolate  spheroid  has  a  low  aspect  ratio  of  l.OOS  with 
=  10.0 .  The  pseudo-spherical  parameters  were  selected  to  facilitate  comparison  to  the  spherical  shell 
which  was  modeled,  tested  and  presented  in  Jones-Oliveira  (Ref.  2).  It  should  be  noted  that  the  parameters 
are  not  identically  those  cited  in  Jones-Oliveira  (Ref.  2);  therefore,  the  results  will  be  different.  The  shell 
and  fluid  properties  differ  slightly;  and  the  shell  thickness  was  determined  to  ensure  neutral  buoyancy. 
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Case  II:  Model  of  an  intermediate  aspect  ratio  geometry.  Specifically,  the  aspect  ratio  is  V2  with  I0V2 
also.  The  parameters  were  selected,  for  comparison  purposes,  to  match  as  closely  as  possible  those  of 
Pauwelussen  (Ref.  7)  and  Nemergut  &  Brand  (Ref.  6). 

Case  III:  Model  of  a  submarine-like  structure.  The  prolate  spheroid  has  a  high  a^)ect  ratio  of  10.0  and 
^  =  l.OOS.  In  particular,  the  geometry  is  similar  to  a  small  submersible  currently  being  developed  at  The 
Charles  Stark  Draper  Laboratory. 

Numerical  results  are  also  presented  for  a  sphere  which  differs  only  slightly  from  the  pseudo- 
^erical  problem.  Its  geometric  parameters  are  those  given  for  Case  I  with  the  following  exceptions:  [1] 
the  dimensional  semifocal  length  is  zero,  o>  [2]  tin  semiminor  axis  length  is  equal  to  the  semimajor 
axis  length,  b  =  a;  and  [3]  the  semimajor  axis  length  is  equal  to  the  limiting  value,  ^/=  a  . 
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Table  2.  Geometric  parameters:  This  table  provides  the  geometric  parameters  defining  the  three 
illustrative  examples.  Case  I  refers  to  the  low  aqiect  ratio  example.  Case  II  refers  to  the 
intermediate  aspect  ratio  example,  and  Case  in  refers  to  the  hi^  a^iect  ratio  example. 


CASE  I 

CASE  II 

CASE  m 

io 

10.0000 

1.4142 

1.005 

f 

0.0323  m 

1.4142  m 

10.9182  m 

a 

0.3239  m 

2.0000  m 

10.9728  m 

h 

0.3222  m 

1.4142  m 

1.0932  m 

h 

0.G147  m 

0.0549  m 

0.0074  m 

h 

4.5469  X  10'^ 

0.0275 

6.7440  X  10"* 

AR 

1.0050 

1.4142 

10.0375 

e 

0.1000 

0.7071 

0.995 

jL. 

\2f^ 

1.7229  X  10'^ 

1.2554  X  lO"* 

3.9  X  10  * 

M 

2.8224 

4.6761 

189.4994 

T 

4,642.9251  t 

749.0  t 

136.5194  t 

/3 

1.6091 

0.2596 

0.0473 

fUnihertz) 

736.073 19£1« 

ll9.20705ft« 

21.72774n. 
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3.2  FREE  VIBRATIONS  OF  PROLATE  SPHEROIDS. 


The  first  step  in  verifying  the  DOE-MACSYMA  coding  of  the  analysis  was  to  compute  the 
free  in  vacuo  vibrations  of  the  three  representative  prolate  spheroids.  Included  are  the  couplings  of  the 
coefficients  associated  with  the  first  eight  Legendre  polynomials.  The  spherical  results  are  presented  for 
comparison  purposes.  The  nondimensionalized  results  are  tabulated  in  Table  3  on  page  12  and  the 
dimensionalized  frequencies  are  tabulated  in  Table  4  on  page  14. 

It  is  noted  that  the  in  vacuo  natural  frequencies  of  the  pseudo-spherical  isolate  spheroid,  i.e.. 
Case  I,  for  which  4  =  10.0  are  within  1%  of  those  computed  for  a  qihere  using  the  same  parameters. 
However,  one  would  expect  the  frequencies  of  the  sphere  to  be  slightly  higher  than  those  of  the  prolate 
spheroid.  The  discrepancies  may  be  explained  by  the  truncation  of  the  infmite  series  lepresentaiions  of  the 
prolate  ^eroidal  angular  functions  expressed  in  terms  ttf  Legendre  polynomials  and/or  by  numerical  errors 
and/or  by  omission  of  the  expansion  coefficients. 

Agreement  of  the  first  few  modes  of  Case  n,  for  which  ^  =  V2,  is  remarkable.  Nemergut  & 
Brand  (Ref.  6)  computed  an  identical  frequency  for  the  lower  branch  of  mode  2.  In  their  paper,  they  compare 
their  own  results  with  those  of  Shiraishi  &  DiMaggio  (Ref.  8)  and  Silbiger  &  DiMaggio  (Ref.  9). 

The  results  for  Case  in  for  which  ^  s  l.OOS,  have  to  be  considered  somewhat  suspect  The  reasons  for 
concern  are:  [1]  the  lowest  odd  frequency  is  not  close  enough  to  zero;  and  [2]  the  branches  cross  each  other. 
Clearly,  the  rigid  body  translation,  for  which  m  =  1,  should  have  a  zero  root  Normalization  by  the  lowest 
"ncMizero"  frequency  reveals  cause  for  ctmcem.  Regarding  the  crossing  of  the  branches,  it  is  not  known 
whether  or  not  these  results  are  correct  Several  options  for  further  investigation  are:  [1]  increase  the 
precision;  [2]  increase  the  number  of  modes  that  are  allowed  to  couple;  and/or  [3]  include  the  expansion 
coefficients  in  the  analysis.^  None  of  these  suggestions  is  numerically  trivial;  however,  the  matter  will  be 
addressed  further. 


^  As  yet  the  expansion  coefficients  in  the  numerics  are  all  taken  to  be  one.  The  assumption  is  that  this 
dependence  will  ultimately  be  absorbed  into  the  time-dep»ident  generalized  displacements  as  per  (Ref.  4, 
Chapter  2.6.3). 
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It  should  be  noted  that  thoe  are,  in  fact,  two  branches;  i.e.,  two  natural  frequencies  associated 
with  each  mode  number.  Please  refer  to  Figure  2  on  page  13,  which  plots  the  dimensional  frequencies  as  a 
function  of  mode  number  for  Cases  I  &  n.  This  phenomenon  is  an  inherent  characteristic  of  axisymmetric 
^ell  structures.  The  simultaneous  solution  of  the  homogeneous  system  of  partial  differential  equations, 
which  governs  the  normal  and  tangential  deflections  of  shells  of  revolution,  reduces  to  the  solution  of  a 
quartic  equation;  specifically,  it  is  bi-quadradc.  Except  fw  the  pure  breathing  mode  m=  0,  there  result  two 
sets  of  nontrivial  roots  to  the  frequency  determinant  resulting  in  a  bifurcation  into  a  lower  branch  and  an 
uiqier  branch  (Ref.  S,  Chapter  8).  The  m  =  0  is  a  special  case;  it  has  only  a  single  complex<onjugate  pair 
of  roots  because  Uo  «  0. 

The  lower  branch  appears  to  contain  the  ~  parameter  and  is,  therefore,  sensitive  to  bending 

12'* 

stiffness.  On  the  other  hand,  the  upper  branch  is  reladvely  insensiuve  to  bending  suffness.  If  one  wer^o 

model  only  the  extensional  vibradons  of  the  prolate  spheroidal  shell,  which  would  imply  that  the  — ^ 

12^* 

terms  were  omitted,  then  the  frequencies  on  the  lower  branch  would  be  too  low  but  the  frequencies  on  the 
upper  branch  would  not  be  affected  appreciably. 

Finally,  all  of  the  even  (odd)  normal  displacements  are  coupled  to  all  of  the  even  (odd) 
tangential  displacements.  Therefwe,  it  is  not  possible,  a  priori,  to  discern  which  frequency  (lower  or  upper 
branch)  is  associated  with  which  displacement  It  is  conjectured  that  the  lower  branch  may  be  dominated  by 
the  tangential  displacements. 
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Table  3.  Nondimensionalized  in  vacuo  natural  frequencies. 


SPHERE 

PROLATE 

SPHEROID 

Case  I 

Case  n 

Case  in 

mode 

T 

^  =  10.0 

§)  =  V2 

5)  =  1.005 

1 

0.0 

0.0 

0.0049 

2.7448 

2 

2.5325 

2.5492 

3.7592 

7.7304 

3 

3.0252 

3.0470 

4.6356 

12.9879 

4 

3.2799 

3.3039 

4.9566 

17.9642 

5 

3.5125 

3.5385 

5.2955 

23.4785 

6 

3.8030 

3.8311 

5.5944 

28.4777 

7 

4.1922 

4.2232 

6.2635 

37.3266 

8 

4.7012 

4.7358 

6.6211 

36.9838 

0 

5.8138 

5.8332 

7.1347 

41.6152 

1 

7.1210 

7.1544 

9.5473 

39.1475 

2 

9.8168 

9.8457 

12.3497 

44.4057 

3 

13.1101 

13.1445 

15.7277 

48. 1'*  ■71 

4 

16J785 

16.6207 

19.5272 

52.5743 

5 

20.1094 

20.1602 

23.5383 

73.3820 

6 

23.6685 

23.7281 

27.6355 

81.2022 

7 

lllATlb 

27.3113 

32.8058 

198.7833 

8 

30.8255 

30.9032 

372074 

218.9494 
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Figure  2.  Dimensional  Frequency  versus  mode  number. 
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Table  4.  Dimensional  in  vacuo  natural  frequencies  (kilohertz). 


SPHERE 

PROLATE 

SPHEROID 

Case  I 

Case  n 

Case  HI 

mode 

r 

§)  =  10.0 

5o  =  V2 

5o  =  1.005 

1 

0.0 

0.0 

5.8412  X  lO'^ 

0.05% 

2 

1.8639 

1.8764 

0.4481 

0.1680 

3 

2.2265 

2.2428 

0.5526 

0.2822 

4 

2.4140 

2.4319 

0.5909 

0.3903 

5 

2.5852 

2.6046 

0.6313 

0.5101 

6 

2.7990 

2.8200 

0.6669 

0.6188 

7 

3.0854 

3.1086 

0.7467 

0.8109 

8 

3.4601 

3.4859 

0.7893 

0.8036 

0 

Aim 

4.2937 

0.8720 

0.9042 

1 

52411 

5.2662 

1.1381 

0.8506 

2 

7.2252 

7.2472 

1.4722 

0.9648 

3 

9.6491 

9.6753 

1.8749 

1.0461 

4 

12.2018 

12.2341 

2.3278 

1.1423 

5 

14.8005 

14.8394 

2.8059 

1.5944 

6 

17.4200 

17.4656 

3.2943 

1.7643 

7 

20.0506 

20.1031 

3.9107 

4.3191 

8 

22.6876 

22.7470 

4.4354 

4.7573 
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3.3  FLUID-LOADED  VIBRATIONS. 


In  order  to  solve  the  even-even/odd-odd  coupled  equations  for  the  fluid-loaded  natural 
frequencies,  it  is  necessary  to  simplify  the  equations.  Here,  the  emphasis  is  on  the  "coupling"  of  the 
coefficients  in  the  differential  equations  resulting  from  the  Legendre  polynomials  representation.  The 
trimming  procedure  is  as  follows: 

1 .  Inspect  the  coefficients  of  each  displacement  ^ipearing  in  the  left-hand  side  of  each 
of  the  modal  equations  given  by  (Ref.  4,  3.S.2a  &  3.S.2b).  Normalize  the 
coefficients  by  the  coefficient  of  the  di^lacement  in  the  n'*  equation.  Let  each 
ncKinalized  coefficient  which  is  less  than  0.0001  be  set  to  0.0.  Allow  an  additional 
order  of  magnitude  for  each  power  of  s  (since  the  observed  range  of  s  satisfles 
li  I  <  20). 

2.  Note  that  the  coupling  terms  at  the  fluid-solid  interface  tqqiearing  on  the  right-hand 
side  of  equation  (Ref.  4,  3.S.2a)  (see  also  Ref.  4,  3.3.5.4)  have  a  factor  of  bni  x 
btim-  ^3ch  of  these  coefficients  drops  off  rapidly  as  the  mode  numbers  increase,  so 
their  products  drop  off  even  faster.  Practically,  this  means  that  the  inner  summation 
in  the  numerator  of  (Ref.  4, 33.S.4)  is  not  infinite;  in  fact,  in  most  cases,  only  the 
diagonal  elements  survive  the  trimming  procedure.  In  other  words,  only  the  {bnnif^ 
terms  usually  survive. 

3.  When  necessary  and  appropriate,  back-substitute  the  solutions  for  the  lower  modes 
into  the  equations  for  the  higho-  modes  to  solve  for  the  higher  modes. 

These  simplifications  are  necessary  in  order  to  compute  the  solution  effectively.  Otherwise, 
the  order  of  the  polynomials  which  must  be  solved  exceeds  the  capabilities  of  the  root  solvers.  For 
example,  consider  Case  1  for  which  ^  =  10.  Without  these  simplifications,  the  order  of  the  polynomial 
for  solving  the  coupling  of  wq.  U2.  H4  >  W4  is  64  and  the  DOE-MACSYMA  ALLROOTS  solver 
was  only  able  to  fund  14  of  the  roots.  The  polynomials  become  ill-conditioned. 
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It  was  verified  that  the  roots  are  quite  stable  against  this  trimming  procedure.  For  the  smaller 
problem,  which  models  the  the  most  significant  cross-coupling,  specifically  that  associated  with  the  lowest 
modes(wni  U2  and  W2),  the  roots  moved  by  less  than  1%  with  these  simplifications.  This  smaller  problem 
could  be  solved  exactly  and  the  results  were  used  for  the  verification. 

The  characteristic  roots  are  presented  in  the  fcnn 

s  -  -  a  ±Qi 

where  s  =  Laplace  variable 
-  <r  =  tenuation 

A  =  nondmensionalized  damped  natural  fie(|iency. 

The  first  few  roots  fOT  Case  I,  for  which  |a  =  10.0,  are 

m  =  0:  Wo,  «o  «  0 

j  =  - 1.940 
s  =  -  1.366  ±  5.400i 

m  =  1;  wi.ui 

s=  0 

j  =  .  1.526  ±0.918i 
j  =  - 0.887+ 6.8301 

m  =  2:  W2,  U2 

j  =  -  1.097 
j  =  -  1.940 
s  =  -  3.052 
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s»-1.9101±2.455i 
j  =  -  0.1306  ±  1.785i 
j  =  - 2.692  ±  1.779/ 
s  =  -  0.403  ±  9.701i 


3.4  SHOCK-LOADED  VIBRATIONS. 


The  characteristic  roots  for  the  shock-loaded  problem  are  identical  to  the  roots  given  for  the 
fluid-loaded  homogeneous  problem  with  the  addition  of  a  root  associated  with  the  shock.  Specifically,  there 
is  an  additional  teal  root  s  =  -1.609. 

The  nondimensionalized  normal  and  tangential  displacements  the  shock-loaded  submerged 
shell  are  plotted  as  a  function  of  nondimensionalized  time.  The  dimensional  deflections  may  be  detennined 
from  the  nondimensionalized  values  by  multiplying  by  the  factors  fm  the  given  geometry:  Case  I  for 
which  ^  =  10.0.  0.323:  Case  II  for  which  ^  =  V2,  %(J=  2.0;  and  Case  III  for  which 

^  =  1.005,  ^7=  10.97 

Similar  conversions  may  be  used  to  obtain  the  physical  dimensional  velocities  and 
accelerations.  Refer  to  Tatde  5  on  page  18  as  necessary  to  relate  nondimensionalized  time  x  to  real  time  t. 
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Table  S.  Nondimensionalized  time  versus  real  time  (ms). 


T 

CaMl 

CaMU 

CaMin 

i  S  10.0 

II 

0 

0.0 

OJ) 

0.0 

1 

0.21S 

1.33S 

7.325 

2 

0.431 

X670 

14.650 

3 

0.646 

4.00S 

21.975 

4 

0.S62 

S.34I 

29.300 

s 

1.0T7 

6.676 

36.625 

10 

X1S4 

13.3S1 

73J50 

IS 

3.231 

20.027 

109475 

20 

4.30S 

26.702 

146.500 

2i 

S.3SS 

33.378 

183.124 

30 

6.461 

40J)S3 

219.749 

3S 

7.538 

46.729 

256.374 

40 

S.61S 

53.405 

29X999 

4S 

9.692 

60.080 

329.624 

SO 

ia769 

66.7S6 

366.248 

60 

1X923 

10.107 

439.498 

70 

13.077 

93.4S8 

51X748 

SO 

17.231 

106.809 

585.997 

90 

19.384 

120.160 

659.247 

100 

21.538 

133.511 

73X497 

no 

23.692 

146.863 

805.763 

i» 

2S.S46 

160.214 

878.996 

130 

28.000 

173.S6S 

95X246 

140 

30.1  S3 

186.916 

1025.49S 

ISO 

32.307 

200.267 

1098.745 
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3.4.1  Case  I:  Low  Aspect  Ratio. 


Figure  3  on  page  21  through  Fgure  16  on  page  34  are  DOE-MACSYMA-generated  plots  of 
the  normal  and  tangential  generalized  displacements  (i.e..  the  time-dependent  amplitudes),  normal  and 
tangential  physical  displacements,  velocities  and  accelerations  for  the  low  aspect  ratio  geometry  in  response 
to  the  qiecified  underwater  explosion. 

The  following  observations  were  made. 

•  The  rigid-body  generalized  displacements  (modal  time-dqiendem  amplitudes),  w  i . 
and  Ui,  indicate  a  finite  displacement  of  the  shell.  This  may  be  explained  as 
follows.  Recall  that  the  generalized  forces  are  the  superposition  of  the  incident, 
scattered  and  radiated  pressure  waves  in  the  fluid.  The  incident  shock  wave  imparts  a 
force  onto  the  prolate  sphermdal  shell  and  continues  past  the  shell.  Scattered  energy 
is  reflected  off  the  shell  back  into  the  fluid.  Radiated  energy  is  also  put  back  into 
the  fluid  due  to  the  motion  of  the  shell  through  the  fluid,  which  can  also  te-Ioad  the 
structure.  The  loss  of  energy  to  the  fluid,  once  the  exponentially  decaying  incident 
shock  wave  has  passed  the  prolate  spheroid  (r  >  2),  damps  the  motion  of  the  shell 
as  it  moves  through  the  fluid  resulting  in  the  finite  displacement 

•  It  is  also  noted  that  max(w2)»  max  (wo)  This  means  that  W2  is  strongly 
coigtled  to  wo  but  not  vice  versa.  The  equations  coupling  wo,  ui.  and  W2  reveal  that 
wo  has  both  of  the  toots  fw  m  =  0  and  the  fourth  root  for  m  =  2,  which  points  out  a 
strong  coupling  to  ui  .  This  can  be  deduced  because  this  fourth  m  =  2  root  is  the 
extra  toot  associated  with  U2,  which  is  predicted  by  the  analysis  and  comparison  to 
the  spherical  problem. 

•  max  (W2)»  max  (wo),  max(w4),  -  andmax(w3)»  max(w5),  max(w7), 
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Ua  drops  off  rapidly  as  n  increases. 


*  As  the  mode  number  increases,  the  couplings  become  negligible. 

*  As  time  progresses,  the  displacements  can  be  seen  to  damp  out  and  then  to  re-excite 
as  the  fluid  and  sheU  motions  exchange  their  energy  through  the  kinematic  boundary 
conditions  at  the  fluid-solid  interface. 

*  The  number  of  modes  required  to  modd  the  reqxxise  adequately  was  determined  by 
examining  the  accelerations,  rather  than  by  considering  just  the  relative  magnitudes 
of  the  di^lacement  coefficients.  It  was  determined  that  the  acceleration  contributed 
by  the  7‘^  mode  was  less  than  10%  of  the  total  acceleration  contributed  by  modes 
m  =  0  -  6.  A  better  criterion  could  be  established  by  examining  the  conservation  of 
energy  of  the  system.  Acceleration  information  is  of  keen  interest  to  experimental 
engineers  because  these  predictions  can  be  verified  directly  using  accekrometers. 

The  shell  strains,  as  defined  by  (Ref.  4, 3.1.13),  and  strain  rates  are  presented  in  Figure  17  on 
page  3S  through  Figure  20  on  page  38.  Again,  the  strain  information  is  of  interest  to  experimental 
engineers  because  strain  gauge  readings  are  easily  obtained.  Comparisons  of  predicted  values  versus 
empirical  data  may  be  accomplished  using  the  square  root  sum  of  the  squares. 
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Legend 


Figure  3. 


-  Wl  W7 

•  W3 


w« 


Modal  time-dependent  nonnal  amplitudes  for  m  =  0-7:  With  the  exception  of  component 
m=l,  the  generalized  displacements  return  to  their  original  equilibrium  position.  However, 
the  m  s  1  contribution  results  in  a  finite  displacement  of  the  shell. 
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Figure  4. 


Transient  normal  deflection  coefficients  fOT  even  m  =  0  -  6;  Note  that  the  m  =  2  mode  has 
the  most  significant  amplitude.  This  explains  its  strong  coupling  to  the  lower  even  modes. 
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Figure  5. 


Transient  normal  deflection  coefficients  for  odd  m  =  1  -  7:  The  m  =  1  mode  finite 
displacement  reflects  the  transfer  of  energy  from  the  shell  to  the  fluid  in  the  form  of 
radiation. 
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Figure  6.  Early-time  normal  displacements  at  five  locations  on  the  shell:  The  positions  vary  from  the 
nose  to  the  back  side  relative  to  the  shock. 
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Figure  7.  Longer-time  normal  di^lacements  at  five  locations  on  the  shell:  Note  the  oscillation  of  the 
displacement  amplitudes — first  dying  out  and  then  re-exciting. 
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Legend 


Figure  8. 
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Transient  normal  velocities  at  five  locations  on  the  shell:  Additional  data  points  would 
smooth  out  the  plots. 
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Legend 
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Transient  normal  accelerations  at  five  locations  on  the  shell:  Additional  points  would 
smooth  out  the  plots. 
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Figure  10.  Modal  time-dependent  tangential  amplitudes  for  m  s  i  -  7;  With  the  exception  of 
component  ms  1,  the  generalized  displacements  return  to  their  original  equilibriu  n 
position.  However,  the  m  =  1  contribution  results  in  a  finite  displacement  of  the  shell. 
Note  that  m  =  0  is  zero  by  definition. 
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Figure  11. 


Transient  tangential  deflection  coefficients  for  even  m  =  2  -  6:  Note  that  the  m  =  2  mode 
has  the  most  significant  amplitude.  This  explains  its  strong  coupling  to  the  lower  even 
modes. 
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Figure  12. 


Transient  tangential  deflection  coefficients  for  odd  m  =  1  •  7:  The  m  =  1  mode  finite 
displacement  reflects  the  transfer  of  energy  from  the  shell  to  the  fluid  in  the  form  of 
radiation. 
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Figure  13. 


Eurly-tiine  tangential  displacements  at  three  locations:  The  tangential  displacements  at  the 
two  vertices  are  zero. 
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Figure  14.  Longer-time  tangential  displacements  at  Hve  locations;  Note  the  oscillation  of  the 
displacement  amplitudes — first  dying  out  and  then  re-exciting,  albeit  to  a  lesser  extent  than 
was  the  case  for  the  normal  displacements. 
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iMgure  IS.  Transient  tangential  velocities  at  three  locations  on  the  shell;  Additional  data  points  would 
smooth  out  the  plots. 
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Transient  tangential  accelerations  at  three  locations:  Additional  data  points  would  smooth 
out  the  plots,  and  perhaps  slightly  increase  the  maximum  values. 
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3.4.2  Case  11  and  Case  III. 


Reasonable  numerical  results  could  not  be  obtained  for  these  cases.  We  believe  that  this  did 
not  result  from  a  fundamental  flaw  in  the  theory.  Rather,  it  is  because  the  series  solutions  do  not  converge 
as  fast  as  the  a.'qwct  ration  is  increased.  In  addition,  cross-coupling  terms  which  were  negligible  in  the  first 
case  are  more  important  in  these  cases.  If  the  model  size  that  DOE-MACSYMA  could  handle  were 
increased  and  more  computer  power  were  used,  solution  should  be  possible.  This  is  not  currently  practical. 
However,  as  computer  power  expands  in  the  future,  it  should  become  more  reasonable. 
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SECTION  4 


CONCLUSIONS 


In  order  to  solve  the  prolate  spheroidal  fluid-solid  interaction  problem,  it  was  necessary  to 
balance  classical  analysis  and  numerical  code  develoinnenL  Separation  of  the  prolate  spheroidal  wave 
equation  reveals  the  prolate  spheroidal  angular  and  radial  functions  to  be  the  exact  solutions  to  the  prolate 
qjhermdal  acoustic  wave  equation,  as  expected.  The  prolate  q>heioidal  angular  functions  are  also  determined 
to  be  eigen  solutions  for  the  prolate  spheroidal  shell  equilibrium  equations.  The  shell  displacement  and 
fluid  pressure  fields  are  exjnessed  in  terms  of  an  infinite  series  expansion  using  the  complete  set  of  prolate 
spheroidal  angular  functions  and  are  therefore  exact  The  resulting  expressions  are  very  complex  (Ref.  4, 
Chapter  S).  However.  DOE-MACSYMA  was  very  helpful  in  dealing  with  the  complex  algebraic 
manipulations. 

The  final  result  was  an  exact  solution  for  the  transient  response  of  both  the  prolate  spheroidal 
shell  and  the  surrounding  fluid  medium.  However,  the  results  are  represented  in  terms  of  infinite  series, 
embedded  in  infinite  series,  embedded  in  infinite  series,  etc.  Computationally,  such  a  representation  is  very 
difficult  to  deal  with  (Ref.  4,  Chapter  3). 

A  novel  trimming  procedure,  which  was  guided  by  the  functional  insight  gained  in  the 
analysis,  was  required  to  solve  fw  numerical  answers.  This  trimming  procedure  was  sufficient  to  obtain 
shock-loaded  results  for  the  nearly-spherical  prolate  spheroidal  geometry  (Ref.  4,  Chapter  S).  Fot  this  case, 
the  results  were  compared  with  the  results  previously  generated  for  the  perfect  sphere  and  the  agreement  was 
excellent 


Unfortunately,  this  procedure  was  not  sufficient  to  obtain  valid  solutions  for  the  higher 
eccentricity  geometries.  One  of  the  primary  problems  for  the  eccentric  geometries  is  modal  coupling. 
Unfortunately,  the  fact  that  the  resulting  series  solutions  are  separately  solutions  to  the  wave  equations  and 
the  shell  equations  does  not  mean  that  they  are  normal  eigenfunctions  which  would  uncouple  the  equations. 
Theoretically,  a  set  of  normal  eigenfunctions  exist  which  would  uncouple  the  equations.  However,  it  is  not 
known  how  to  find  such  a  set  of  functions. 

In  principal,  the  equations  can  be  evaluated  in  spite  of  the  cross-coupling.  However,  it  was 
determined  that  the  higher  the  eccentricity,  the  slower  the  power  series  converged  and  the  more  significant 
the  coupling  of  the  modes.  For  practical  solution  with  currently  reasonable  computer  resources,  the  cross¬ 
coupling  terms  must  be  dropped  and  the  series  truncated  at  a  reasonable  number  of  terms.  Unfortunately, 
this  results  in  incorrect  answers.  In  the  future,  increased  computer  power  may  make  this 
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solution  useful  in  an  engineering  sense.  However,  at  this  time  the  method  can  not  be  used  to  solve 
practical  problems. 

Since  reasonable  numerical  solutions  could  not  be  obtained  for  high  eccentricity  geometries 
using  the  new  closed-form  theory,  we  were  unable  to  compare  the  new  theory  with  experimental  or  fmite 
element  results. 
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